We investigate signatures of helium diffusion in late F-type 
stars by astero seismology. Stellar models were computed with 
different physical inputs (with or without element diffusion) 
and iterated in order to fit close-by evolutionary tracks for each 
mass. The theoretical oscillation frequencies were computed 
and compared for pairs of models along the tracks. Various as- 
teroseismic tests (large separations, small separations, second 
differences) were used and studied for the comparisons. The 
results show that element diffusion leads to changes in the fre- 
quencies for masses larger than 1.2 M . In particular the he- 
lium gradient below the convective zone should be detectable 
through the second differences. 

Key words, asteroseismology; diffusion; stars: oscillations 
f-*^ (including pulsations); stars: interiors 
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Abstract. Element diffusion is expected to occur in all kinds of stars : according to the relative effect of gravitation and 
radiative acceleration, they can fall or be pushed up in the atmospheres. Helium sinks in all cases, thereby creating a gradient 
at the bottom of the convective zones. This can have important consequences for the sound velocity, as has been proved in the 
sun with helioseismology. 



1. Introduction 

Among the different physical processes that occur in stellar in- 
teriors, element diffusion plays an important role. Under the 
oppositing effects of gravitation and thermal settling on the one 
hand, and radiative acceleration on the other hand, the vari- 
ous elements un dergo relative sep aration. This process, which 
was described by Michaud ( 1970) to account for the abundance 
anomalies observed in chemically peculiar A stars, is now rec- 
ogniz ed to occur in all kinds o f stars. 

IVauclair & Theadol J2004I described in terms of the so- 
called "second differences" an asteroseismic signature of he- 
lium diff usion in m ain-sequence A stars (1.6 M and 2.0 M Q ). 
iTheado et alJ Mi presented a first discussion of asteroseis- 
mic signatures of diffusion in stars between 1.1 and 1.3 M Q . 
In the present paper we show evidence of differences in the 
internal structure of stars between 1.1 and 1.45 M due to 
element diffusion. We compute evolutionary tracks of stellar 
models with similar external parameters (luminosities, effec- 
tive temperatures and chemical compositions), with or without 
helium diffusion. We then study asteroseismic tests, particu- 
larly the signatures of helium gradients below the outer con- 
vective zones. 

In Section 2, we present the model computations and cali- 
brations. The asteroseismic tests of internal structure between 
models with and without diffusion are discussed in Section 3. 
The helium gradients and their consequences for the second 
differences are studied in Section 4. Section 5 gives the discus- 
sion and conclusion. 

2. Calibration and characteristics of the models 

We computed models with the Toulo use-Geneva Evolutionary 
Code (TGEC) jRichard et alJll996l) . for masses from 1.1 to 
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1.45 M , respectively noted Ml.l to Ml. 45 (see Tabled). The 
physical inputs and parameters are : 

- Equation of state : MHD (Dapp erJ 19921 

- Opacities : OPAL (Iglesias & Rogers 1996) completed by 
the lAlexander & Fergussonl ( ll994l) low temperature opaci- 
ties, 

- Nuclear rea ction rates : NACRE compilation of nuclear re- 
action rates ( Angul oet al J 19991) with the Bahcall screening 
routine, 

- Convection treatm ent : mixing-length theory 
jBohm-Vitensell958l) . 

- Diffusion : diffusion coefficients computed as in 

IPaauetteet alJJl986h . 

For a given stellar mass, we computed two series of mod- 
els : with and without element diffusion. No mixing process or 
mass loss is taken into account. The series were calibrated to 
obtain evolutionary tracks very close in the HR diagram, such 
that, for the same luminosity, their temperature differences are 
never larger than 60 K. This calibration is achieved by adjust- 
ing two free parameters of the stellar evolution code : the initial 
helium abundance Yo and the mixing length parameter a. We 
used for all the homogeneous models the same values of the 
free parameters a and Yo. For each mass, models including 
diffusion are then calibrated to have evolutionary tracks very 
close to the homogeneous ones. The calibration parameters are 
presented in Tabled 

Figure ^ displays the computed evolutionary tracks. For 
each mass the tracks obtained with and without diffusion lie 
very close to each other. Furthermore, up to 1.35 M the age 
evolution is similar in the two cases, so that models with the 
same age lie close together in the HR diagram. Regarding the 
uncertainties, we can consider that these models are observa- 
tionaly identical, i.e. with the same effective temperatures, lu- 
minosities and chemical compositions. This is no longer the 
case for larger masses : along track Ml. 45, of 1.45 M , the 
ages for models with diffusion are quite different to those with- 
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Models 


Mass (Mo) a Y 


Ml.l-hom 
Ml.l-dif 


1.1 1.75 0.268 
1.1 1.84 0.267 


M1.2-hom 
M1.2-dif 


1.2 1.75 0.268 
1.2 1.88 0.268 


M1.3-hom 
M1.3-dif 


1.3 1.75 0.268 
1.3 1.99 0.268 


M1.35-hom 
M1.35-dif 


1.35 1.75 0.268 
1.35 1.99 0.268 


M1.45-hom 
M1.45-dif 


1.45 1.75 0.268 
1.45 1.88 0.278 




3.74 



Fig. 1. Evolutionary tracks of the 1.1 to 1.45 M models. 
Dashed lines: homogeneous models, solid lines: models includ- 
ing microscopic diffusion 



out diffusion. Indeed, for masses greater than 1.3 M , more 
than half of the energy production is dominated by the CNO 
cycle, which is more sensitive to the central temperature than 
the pp chains. Since models with diffusion have a larger initial 
helium fraction than models without diffusion, the former are 
more evolved for the same age. 

Table |2] gives the characteristics of the models with and 
without diffusion for each mass and different ages. The cali- 
brated models do not have strictly identical values of the sur- 
face parameters but they are very close, except for the case 
Ml. 45 for which models with similar parameters have differ- 
ent ages. 



3. Tests of internal structure 

The oscillation frequencies of our models have been computed 
with an updated v ersion of the adiabatic code described in 
iBrassard et al.l dl992l) . for values of the azimuthal degree I be- 
tween and 3 and different values of the radial order n. We then 
studied in detail different combinations of these frequencies 
which could lead to observational tests of the internal structure 
and chemical composition of the stars. 

The "large separations" represent the frequency differences 
between two modes of the same degree / and successive radial 
orders n: 

Av„,; = V n ,l ~ Vn-\,l (1) 

For p-modes in the asymptotic theory (n » I), the large sep- 
arations are nearly constant. The so-called "echelle diagrams" 
present the frequencies in ordinates, and the same frequencies 
modulo the average large separation in absissae. The asymp- 
totic theory predicts a vertical line for each degree. 

The differences between the lines / = and / = 2 on the 
one hand, / = 1 and I — 3 on the other hand, for each case, 
are the so-called "smal l separations" , whic h we normalize as 
suggested by Roxburgh & Vorontsovl 11200 1 ) : 

Dv u+2 = 2 ^ (v„j - v„_ u+2 ) (2) 

The deviations of the small separations to zero show departure 
from the asymptotic theory. The small separations are mainly 
sensitive to the stellar core. 

In Figure |2] we present as an example the echelle diagram 
for models Ml.l to Ml. 35 at different ages. For Ml. 45, we 
compare echelle diagrams for models at the same location on 
the HR diagram, which correspond to different ages for the 
models with and without diffusion. This is presented in Figure 
[3] The small separations of the same models are presented in 
Figure@]for models Ml.l to Ml .35 and in Figure|5]for models 
M1.45. 

We can see on these figures that the differences between the 
models with and without diffusion become more important dur- 
ing the evolution when the mass of the star increases. Up to 1.3 
M Q , the tracks in the large and small separations lie very close 
to each other whereas for higher masses significant deviations 
do appear. These effects are related to the internal differences 
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Table 2. Physical characteristics of the models; L/L : luminosity in solar unity; |AL|/L : relative difference of luminosity between 
models with and without diffusion; T e ff : effective temperature; |AT|/T : relative difference of effective temperature between 
models with and without diffusion; X s and Y s : surface hydrogen and helium mass fractions 



Models 


Ages (Gyr) 


L/L 


IALI/L 


Teff (K) 


|AT|/T 


X s 


Y s 


Ml.l-hom 


1.5 


1.26637 




6034.90 




0.7130 


0.2680 




2.2 


1.35947 




6056.06 




0.7130 


0.2680 




4.5 


1.71937 




6092.56 




0.7130 


0.2680 




5.2 


1.86054 




6076.73 




0.7130 


0.2680 


Ml.l-dif 


1.5 


1.26488 


0.12% 


6049.23 


0.24% 


0.7298 


0.2517 




2.2 


1.36013 


0.05% 


6062.48 


0.11% 


0.7379 


0.2439 




4.5 


1.72119 


0.11% 


6071.28 


0.35% 


0.7604 


0.2222 




5.2 


1.86513 


0.25% 


6048.67 


0.46% 


0.7667 


0.2162 


M1.2-hom 


1 


1.89069 




6301.30 




0.7130 


0.2680 




1.5 


2.00586 




6312.19 




0.7130 


0.2680 




2.6 


2.30632 




6316.55 




0.7130 


0.2680 




3.3 


2.48519 




6271.48 




0.7130 


0.2680 


M1.2-dif 


1 


1.89099 


0.016% 


6332.72 


0.50% 


0.7335 


0.2482 




1.5 


2.00812 


0.11% 


6335.34 


0.37% 


0.7438 


0.2382 




2.6 


2.30420 


0.09% 


6315.97 


0.01% 


0.7661 


0.2168 




3.3 


2.47611 


0.37% 


6250.72 


0.33% 


0.7745 


0.2087 


M1.3-hom 


0.5 


2.65461 




6551.19 




0.7130 


0.2680 




1.2 


2.94286 




6556.77 




0.7130 


0.2680 




2.2 


3.33212 




6479.03 




0.7130 


0.2680 




2.6 


3.44310 




6406.05 




0.7130 


0.2680 


M1.3-dif 


0.5 


2.66146 


0.26% 


6615.92 


0.99% 


0.7390 


0.2430 




1.2 


2.95386 


0.37% 


6592.50 


0.54% 


0.7756 


0.2077 




2.2 


3.33165 


0.014% 


6483.36 


0.067% 


0.8033 


0.1810 




2.6 


3.43819 


0.14% 


6411.80 


0.09% 


0.7936 


0.1901 


M1.35-hom 


0.5 


3.18750 




6681.75 




0.7130 


0.2680 




1 


3.42768 




6678.21 




0.7130 


0.2680 




1.8 


3.84282 




6588.55 




0.7130 


0.2680 




2.2 


3.98263 




6495.91 




0.7130 


0.2680 


M1.35-dif 


0.5 


3.19352 


0.19% 


6726.98 


0.68% 


0.7637 


0.2191 




1 


3.43558 


0.23% 


6698.23 


0.30% 


0.8075 


0.1771 




1.8 


3.84034 


0.06% 


6570.07 


0.28% 


0.8429 


0.1428 




2.2 


3.97924 


0.08% 


6481.87 


0.22% 


0.8214 


0.1631 


M1.45-hom 


1 


4.82217 




6906.53 




0.7130 


0.2680 




1.2 


4.98999 




6842.74 




0.7130 


0.2680 




1.6 


5.23299 




6709.96 




0.7130 


0.2680 




1.8 


5.33777 




6607.09 




0.7130 


0.2680 


M1.45-dif 


1 


5.13913 


6.57% 


6769.24 


1.99% 


0.9948 


0.0010 




1.2 


5.32206 


6.65% 


6676.52 


2.43% 


0.9943 


0.0014 




1.6 


5.58805 


6.35% 


6469.64 


3.71% 


0.9823 


0.0104 




1.8 


5.71755 


7.11% 


6347.46 


3.93% 


0.9660 


0.0250 



in chemical composition and to the different extension of the 
convective cores. 



According to the results of iBazot et all J2005I) . which 
presents the seismic analysis of the star jj. Arae, observed with 
the spectrograph HARPS, the uncertainties on the small sep- 
arations are about 0.37 yuHz. In the cases of 1.1 and 1.2 M Q , 
the differences between models with and without helium dif- 
fusion are too small to be observable. For higher masses, this 
difference could be detectable. 



4. Tests of helium gradients 

4.1. The helium gradients 

The largest differences between the models with and without 
diffusion lie just below the convective zone where helium drifts 
inward due to the effect of gravitation. Figures [6] to [8] display 
the helium profiles and the gradients of the sound velocity in 
models with (solid lines) and without (dashed lines) diffusion 
for Ml. 2, Ml. 3 and Ml. 45, as a function of the acoustic depth 
(time needed for the waves to travel from the surface down to 
the considered layer), at different ages. Models Ml.l behave in 
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Fig. 2. Echelle diagram for models Ml.l to Ml. 35 at different 
ages. Crosses : I - 0, stars : I = 1, triangles : / = 2, diamonds 
: I = 3. The points connected by lines are for models with 
diffusion. 



Fig. 4. Small separations for models Ml . 1 to Ml .35 at different 
ages. The lines are for models with diffusion (solid lines : Dv%2, 
dotted lines : Dv\^) and the symbols are for models without 
diffusion (triangles : Dvn,2, diamonds : Dv\^). 




Fig. 3. Echelle diagram for models Ml. 45 at 1.2 Gyrs (with 
diffusion) and 1 .6 Gyrs (without diffusion) on the left, and at 
0.7 Gyrs (with diffusion) and 1 .2 Gyrs (without diffusion) on 
the right. Crosses : I = 0, stars : I = \, triangles : / = 2, dia- 
monds : I = 3. The points connected by lines are for models 
with diffusion. 



Fig. 5. Small separations for models Ml. 45 at 1.2 Gyrs (with 
diffusion) and 1.6 Gyrs (without diffusion) on the left, and at 
0.7 Gyrs (with diffusion) and 1 .2 Gyrs (without diffusion) on 
the right. The lines are for models with diffusion (solid lines 
: Dvn,2, dotted lines : Dv\^) and the symbols are for models 
without diffusion (triangles : Dvq,2, diamonds : DV13). 
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a way similar to models Ml. 2 and models Ml. 35 are similar to 
Ml. 3 : they are available in the electronic version only. In the 
helium profiles, helium gradients due to diffusion are clearly 
seen. In models Ml.l and Ml. 2, this helium gradient remains 
smooth throughout the evolution in the main sequence, while 
in the case of models Ml. 3 to Ml. 45, it becomes sharper with 
time. For example, in models Ml. 3, the bottom of the convec- 
tive zone deepens in the stellar interior around 2.5 Gyrs, which 
explains the discontinuity on the helium gradient and the in- 
crease of the helium mass fraction in Table |2] for the corre- 
sponding model, while in models Ml. 2, the convective zone 
sinks later and slower. Actually, the helium gradients in mod- 
els Ml . 1 and Ml .2 become steeper after the turnoff. In the later 
(Figure|5}, the features around 700 s in the sound speed profile 
are due to the helium ionisation zones while the bottom of the 
convective zone leads to characteristic features at t=1900 s (1 
Gyrs), t=2000 s (1.5 Gyrs), t=2200s (2.6 Gyrs) and t=2500 
s (3.3 Gyrs). We can see in this example that the differences 
between the models with and without diffusion increase with 
time, and that this behaviour increases with increasing stellar 
mass. 

4.2. The second differences 

Asteroseismology provides tools to probe the base of the 
convective zone and the helium gradient. As we have seen, 
in these locations the sound velocity undergoes rapid varia- 
tions, which involve partial reflections of the sound waves. 
These reflections modulate the computed frequencies of 
the stellar spectrum, which i s bette r seen in the so-ca lled 
"second differences" JGo ughl Il9 90t iMonteiro & Thompsonl 
1 998: 1 iRoxburgh & Vorontsovll200ll IVauclair & Theadoll2004l 
The ado et alJ l2005). The second differences are defined as fol- 
lows: 



62V 



Vn-l 



2v„. 



(3) 



In order to identify the different components which modulate 
the oscillations, we have computed the Fourier transform of 
the second differences. The positions of the peaks correspond 
to the modulation periods of the second differences due to the 
partial reflections of the waves, which are equal to twice the 
acoustic dept h of the correspond ing features in the sound ve- 
locity dVauclair & Theadoll2004l) . In Figures [5] to ^2 we can 
then identify for each model the peak due to the helium ion- 
isation zone (the closest to the surface) and those due to the 
bottom of the convective zone, where the helium gradient lies. 
The differences between models with and without diffusion are 
clearly visible. 

In all the cases with diffusion ("Figures l9"l to II It . the peaks 
due to the helium ionisation zones decrease in amplitude during 
evolution, because of the helium concentration decreases in the 
convective zone due to diffusion. 

In the 1 . 1 M and 1 .2 M stars, the diffusion process works 
on a time scale too long to be visible in the helium distribution 
during the main sequence. The amplitude of the peak due to the 
base of the convective zone slightly decreases because of the 
sinking of this zone (see Figure|9j- For higher mass stars with 
diffusion (Figures ITOl and [TO . the peak due to the base of the 
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Fig. 6. Helium profiles (left column) and gradients of the sound 
velocity (right column) in the models with (solid lines) and 
without (dashed lines) element diffusion, for 1 .2 M as a func- 
tion of the acoustic depth at 1, 1.5, 2.6 and 3.3 Gyrs 



convective zone undergoes a strong increase of amplitude when 
the convective zone becomes deeper and steepens the helium 
gradient. In the homogeneous models, there is no increase in 
the amplitude of this peak. This behaviour is clearly due to the 
diffusion process. 



5. Discussion 

The large and small separations represent asteroseismic tests 
of the different internal structure in models with and without 
diffusion. Our computations show that these differences in- 
crease during stellar evolution and that it becomes greater for 
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Fig. 7. Same as Fig.gjfor 1.3 M Q at 0.5, 1.2, 2.2 and 2.6 Gyrs Fig. 8. Same as Fig.|6]for 1.45 M Q at 1, 1.2, 1.6 and 1.8 Gyrs 



increasing mass. This behaviour should be detectable for stellar 
masses larger than 1.2 M Q . 

The behaviour of the second differences of the models are 
clear : due to the effect of diffusion, helium falls below the 
outer convective zone and forms a gradient that increases dur- 
ing the evolution. But, depending on the mass of the star, the 
behaviour of the sound waves is different. For low masses (1.1 
and 1.2 M ), the diffusion process is slow, and the depth of the 
convective zone increases slowly. Hence, the helium gradient 
remains smooth during the main sequence and becomes steep 
only at the turnoff. When the mass of the star is large enough, 
the convective zone deepens more quickly and the helium gra- 
dient steepens more rapidly. There is a strong reflection of the 
pressure waves in the region of the helium gradient, which ex- 
plains the large increase in the amplitude of the peak related to 



the base of the convective zone. Later, helium dilutes inside the 
convective zone, so that the helium concentration, which was 
very low after rapid diffusion, increases again (see Table |2j. 
However, the helium gradient remains steep and the amplitude 
of the peak remains high (see FigureslTUlandim. 

The models presented here have been computed with pure 
gravitational settling. In real stars, the variation in the helium 
abundance and the steepening of the helium gradient resulting 
from diffusion will be smaller, due to the mixing processes be- 
low the convective zone and possible mass loss. We can how- 
ever expect that such a signature of helium diffusion may be de- 
tected in stars slightly more massive than the Sun with CoRoT. 
The accuracy of the instrument is expected to reach 0.1 //Hz 
and should allow detection of peaks in the Fourier transforms 
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Fig. 9. Second differences (left column), Fourier transform of Fig. 10. Same as Fig.|5]for models Ml. 3 of 1.3 M Q at 0.5, 1.2, 

the second differences (right column) of models Ml. 2 of 1.2 2.2 and 2.6 Gyrs 
M Q at 1, 1.5, 2.6 and 3.3 Gyrs, with (solid lines) and without 
(dashed lines) diffusion. 



of the second differences of the oscillation frequencies as a sig- 
nature of helium diffusion. 
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